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Abstract 

We present detailed, low temperature, magnetoresistance and specific heat data of single crystal 
YbNiSi3 measured in magnetic field applied along the easy magnetic axis, H\\b. An initially 
antiferromagnetic ground state changes into a field-induced metamagnetic phase at ~ 16 kOe 
(T — > 0). On further increase of magnetic field, magnetic order is suppressed at ~ 85 kOe. The 
functional behaviors of the resistivity and specific heat are discussed in comparison with those of 
the few other stoichiometric, heavy fermion compounds with established field-induced quantum 
critical point. 

PACS numbers: 72.15.Qm,75.20.Hr,75.30.Kz 
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The "quantum criticality conundrum"- continues to be a focus of attention for many 
theorists and experimentalists.-^^^ In the particular case of magnetic intermetallic com- 
pounds that manifest heavy fermion ground states, magnetic ordering temperatures can be 
tuned by some control parameter to T = 0, bringing the material to a quantum critical 
point (QCP). This tuning is possible because of a delicate balance between the Ruderman 
- Kittel - Kasuya - Yoshida (RKKY) interaction, that defines the magnetic ordering tem- 
perature, and the Kondo effect, that causes screening of the magnetic moments and favors 
non-magnetic ground state.- Recently, in addition to traditionally used pressure and chemi- 
cal substitution, magnetic field has also emerged as a potential control parameter having 
an advantage of being a continuous and common parameter for a number of experimental 
techniques. However, the applicability of a single theoretical description for a QCP reached 
by using different tuning parameters, is still under discussion. 

So far the number of stoichiometric materials exhibiting a field-induced QCP is rather 
small, with only few of them being Yb-based: YbRli2Si2,— YbAgGe,— YbPtln,—. The 
release of entropy associated with the magnetic ordering (which can be used as a rough 
caliper of the sise of the Yb magnetic moment in the ordered state) in the three materials 
ranges from ~ 0.0LRln2 for YbRh 2 Si 2 , to < O.LRln2 for YbAgGe, to ~ 0.6i?ln2 for 
YbPtln. Nevertheless the H — T phase diagrams and functional behavior of a number of 
physical parameters in the vicinity of the field-induced QCP are very similar for all three of 
these compounds. Recently investigated YbNiSis, 13 an orthorhombic, moderately - heavy 
fermion compound with the Neel temperature, T/v = 5.1 K, Sommerfeld coefficient, 7 m 190 
mJ/mol K 2 , and entropy associated with the magnetic ordering ~ 0.6i?ln2, seems to be 
a suitable candidate for further studies of the possible field-induced quantum criticality in 
stoichiometric, Yb-based, heavy fermions. 

Single crystals of YbNiSi3 were grown from Sn flux (see Ref. [l3| for more details). Flux 
residue from the surface of the samples was polished or/and etched out to exclude the effect of 
elemental Sn on measured electrical transport properties. DC magnetization was measured 
up to 55 kOe and down to 1.8 K in a Quantum Design MPMS-5 SQUID magnetometer. Heat 
capacity and standard, 4-probe, AC (/ = 16 Hz) resistivity in zero and applied field up to 
140 kOe were measured in a Quantum Design PPMS-14 instrument using a 3 He refrigerator 
with heat capacity and ACT options respectively. In all measurements an external field was 
applied along easy magnetic axis (if ||6).— Resistivity in applied field was measured in a 
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transverse, H _L I, configuration. 

The low field magnetic susceptibility, \ = M/H, measurements and 2 K magnetization, 
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. The feature associated 



M(H), isotherms are consistent with the data reported in Ref. 
with the antiferromagnetic transition in \ ( or i n d(xT) / dT) 1A shifts to lower temperatures 
with increase of the magnetic field applied (Fig. [T]). The lower field metamagnetic transition 
(Fig. [2]) moves to lower fields with increase of temperature at which the M(H) data were 
taken, at higher temperatures the second feature moves to accessible field range (upper 
left inset to Fig. [2]); this feature corresponds to high field break in M(H) slope reported 
previously^. 

Field- dependent resistivity isotherms are shown in Fig. [3J Two features (marked with 
dashed lines), consistent with those seen in aforementioned magnetization (M(H)) measure- 
ments, are easily recognizable and can be monitored in the H—T domain of the experiments. 
The position of the lower field feature is almost temperature- independent, while the critical 
field associated with the higher field feature decreases with an increase of temperature. 

Representative-temperature dependent resistivity, p(T), curves are shown in FiglU No 
traces of the signal from residual Sn flux are seen in H = data. The resistivity at base 
temperature, T « 400 mK, in zero field is ~ 2pQ cm, consistent with p = cm 
obtained from p = po + AT 2 fit in [kJ. A clear break of the p(T) slope and the lower 
temperature decrease of resistivity below the magnetic ordering temperature, T N , associated 
with a loss of spin-disorder scattering is observed up to 70 kOe. At higher fields no signature 
of magnetic transition can be seen. Moreover, there in no clear signature of traditional form 
of non-Fermi-liquid behavior (linear or close-to-linear temperature dependence of resistivity) 
at fields close to the one at which T N is suppressed to T = (Fig. HJ inset). 

The temperature-dependent heat capacity is shown in Fig. [5j Results in zero field are 
consistent with the published data.— Consistent to aforesaid, the peak in C P (T) moves down 
with an increase of the applied magnetic field, it's height decreases, and it is not recognized 
in the data at 80 kOe and above. Closer examination of H = 70 kOe curve (Fig. [5b) reveals 
that the peak associated with the magnetic long range order is located on a background with 
a broad maximum. It is seen with more clarity if C p /T is plotted as a function of temperature 
(inset to Fig. [5b). This broad maximum moves up in temperature and acquires structure 
with an increase of applied field: the data for H = 120 and 140 kOe can be deconvoluted 
to up to three broad peaks. For H < 50 kOe these broad non-monotonic background is 
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probably obscured by the feature in C P (T) associated with the magnetic transition. The 
plausible explanation for these broad features is that they are Schottky-like contributions 
to the specific heat: the the degeneracy of crystal-electric-field-defined lower energy levels 
of YbNiSi3 (e.g. a quadruplet or two closely spaced doublets) is lifted by the Zeeman term 
of the hamiltonian, and the energy difference between different levels increases further on 
increase of the applied field. This explanation is in in dispute with the assumption of the 
doublet ground state of YbNiSiy^ and calls for additional studies. 

In order to evaluate the change of the electronic specific heat coefficient, 7, heat capacity 
data are re-plotted as C p /T vs. T 2 in Fig. El The small upturn at low temperatures (left 
inset to Fig. E]) that becomes more pronounced at higher fields is probably associated with 
the nuclear Schottky contribution to the specific heat. The Sommerfeld coefficient in the 
specific heat, 7, can be estimated by extrapolating the high temperature linear part of the 
C p /T vs. T 2 curves to T = (dashed line in Fig. ED- As a result of such procedure 7 « 250 
mJ/mol K 2 , practically independent of the applied field can be assessed. This value is close 



to the estimate in Ref. 
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for H = 0. However, this estimate exceeds significantly the 
values of C p /T measured in T — > regime for magnetic fields in which the magnetic order 
is suppressed, pointing out that, at least for 70 kOe < H < 140 kOe, high temperature 
extrapolation of C p /T vs. T 2 is not an unique estimate of 7. This brings back the question 
of the estimate of the entropy associated with the magnetic transition^ and may point out 
that (in lieu of the data on non-magnetic analogue) the non-magnetic contribution, taken as 
extrapolation from above 10 K, was overestimated and the cited value of AS mag ~ 0.6i?ln2 
was underestimated. 

Operationally, one can take the extrapolation to T — > of the lowest temperature 
measurements plotted as C p /T vs. T 2 , see left inset to Fig. El (ignoring the nuclear Schottky 
contribution), as an estimate of 7. The obtained values (right inset to Fig. ED seem to 
follow the trend generally observed in materials close to QCP— iiiii2*I£ii£. However, at least 
in the particular case of YbNiSi3, one cannot discard the possibility that the shifts of the 
Schottky-like peaks to higher temperatures with an increase of the applied field leads to 
a decreasing contribution (since it is taken further from the maximum) to the apparent 
7 = C p /T\t^o, and genuine behavior of the Sommerfeld coefficient is masked by this fastly 
changing contribution. 
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The thermodynamic and transport measurements on YbNiSis in magnetic field applied 
along the 6-axis define the H — T phase diagram (Fig. [7]). This phase diagram is consistent 
(in overlapping H — T domains) with the one recently suggested^ from magnetoresistance 
isotherms for T > 1.85 K. The phase diagram is rather simple, with only two phases in 
the ordered state, and a triple point. The positions of the Schottky-like anomalies are also 
shown, depending on the applied field in a close to linear fashion. A magnetic field of ~ 85 
kOe suppresses the long range order, making the phase diagram look similarly to the cases of 
field-induced QCP recently discussed in several materials.— »ILi2iiiLi& Further understanding 
of the low temperature state of YbNiSis in applied field can be approached by detailed 
analysis of the magnetotransport and thermodynamic data. 

Low temperature resistivity (p(T)) data measured in magnetic field were fitted with 
the equation p(T) = po + AT 13 (as e.g. in Ref. [3) with three fitting parameters: residual 
resistivity, po, pre-factor, A, and exponent, (3. Two sets of fits were performed: from 3 K 
down to base temperature (~ 0.4 K) and from 1.5 K down to base temperature. Both sets of 
fits are giving similar results (Fig. [8]). All three fitting parameters have features associated 
with the critical fields, the lower critical field being essentially temperature-independent 
below 3 K and the upper one only shifting slightly as T decreases from 3 K to 1.5 K (as 
can be seen in A and f3). It is noteworthy that the exponent f3 as calculated from the data 
above ~ 0.4 K, does not fall below (3 = 2, and the changes in f3 as obtained from the 0.4 
K > T > 1.5 K are rather small (Fig. [U upper panel). This is clearly different from the 
case of field-induced QCP in YbRli2Si2,— ^ YbAgGe^ 1 ^ and other materials, although it 
is possible that (3 approaches 1.0 for T < 0.4 K, and in a very narrow field range. 

In summary, thermodynamic and transport measurements on YbNiSis {H\\b) reveal H—T 
phase diagram with two distinct, antiferromagnetic and metamagnetic, ordered phases and 
a triple point. An applied magnetic field of ~ 85 kOe suppresses the magnetic ordering 
temperature. As distinct from several known heavy fermions with field induced QCP, above 
400 niK no non-Fermi-signature in resistivity was observed in the vicinity of the critical 
field for T > 0.4 K, although they may well emerge for T < 0.4 K data. Heat capacity 
measurements suggest that applied magnetic field splits degenerate ground state of YbNiSi 3 
and this field-dependent split of the energy levels manifest itself as Schottky-like anomalies. 
Additionally, possible artifacts in assessing of the Sommerfeld coefficient from heat capacity 
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measurements, that have relevance beyond this particular compound, were disclosed. 
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FIG. 1: Low temperature DC magnetic susceptibility, M(T)/H, measured in YbNiSi3 in different 
applied fields, H\\b. Inset: d(xT)/dT corresponding to the data in the main panel. Gap in the 
data near between 4.3 K and 4.6 K corresponds to the region of unstable temperature control in 
our MPMS-5 magnetometer. 
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FIG. 2: Representative magnetic isotherms, M(H), measured in YbNiSi3 at different temperatures 
for H\\b. Insets: enlarged low- and high-field plots of corresponding dM/dH. Arrows on upper left 
inset mark second, high field, feature, that broadens beyond resolution for T = 5 K. 
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FIG. 3: Low temperature resistivity isotherms, p(H), measured at (from the bottom to the top) 
0.4, 0.76, 1, 1.25, 1.5, 1.75, 2, 2.5, 3, and 4 K. Dashed lines are guide for the eye. 
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FIG. 4: Representative low temperature resistivity curves measured in zero and applied magnetic 
field. Curves are shifted down along the Y-axis by multiplicative of 5 cm for clarity. Numbers 
correspond to applied field in kOe. 
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FIG. 5: Low temperature heat capacity, C P (T), measured in different applied magnetic fields. Inset 
to panel (b): C p /T as a function of temperature. Note that the curve for H = 70 kOe is shown on 
both panels. 
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FIG. 6: Heat capacity in applied field for 70 kOe < H < 140 kOe plotted as C p /T vs. T 2 . Dashed 
line - linear extrapolation from the high temperature data. Insets: left - enlarged low temperature 
part; right - estimate of 7 from extrapolated low temperature linear portion of the C p /T vs. T 2 
graph. 
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FIG. 7: Tentative H — T phase diagram of YbNiSi3 for H\ b. Different symbols correspond to 
the phase lines obtained from different thermodynamic and transport measurements: open circle 
- M(T), filled circle - M(H), open star - C P (T), filled star - p(T), diamond and asterisk - p(H). 
Crosses correspond to the position of broad maxima in C p /T, dashed lines are guides for the eye. 
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FIG. 8: Results of fitting of low temperature resistivity to the form p(T) = po + AT@; units of 
p(T) and po are pfl cm. Circles: fit from 3 K to the base temperature, triangles: fit from 1.5 K to 
the base temperature. Nominal error bars of the fits are shown. Vertical dashed lines mark T — > 
critical fields from the phase diagram above. 
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